home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SGBSL.z / SGBSL
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGBBBBSSSSLLLL((((3333FFFF))))                                                            SSSSGGGGBBBBSSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGBSL   - SGBSL solves the real band system A * X = B  or  TRANS(A) * X =
  10.      B using the factors computed by SBGCO or SGBFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE SGBSL(ABD,LDA,N,ML,MU,IPVT,B,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAABBBBDDDD REAL(LDA, N)
  20.         the output from SBGCO or SGBFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array  ABD .
  24.  
  25.      NNNN INTEGER
  26.         the order of the original matrix.
  27.  
  28.      MMMMLLLL INTEGER
  29.         number of diagonals below the main diagonal.
  30.  
  31.      MMMMUUUU INTEGER
  32.         number of diagonals above the main diagonal.
  33.  
  34.      IIIIPPPPVVVVTTTT INTEGER(N)
  35.         the pivot vector from SBGCO or SGBFA.
  36.  
  37.      BBBB REAL(N)
  38.         the right hand side vector.
  39.  
  40.      JJJJOOOOBBBB INTEGER
  41.         = 0         to solve  A*X = B ,
  42.         = nonzero   to solve  TRANS(A)*X = B , where
  43.         TRANS(A)  is the transpose.  On Return
  44.  
  45.      BBBB the solution vector  X . Error Condition
  46.  
  47.      AAAA division by zero will occur if the input factor contains a zero on the
  48.      diagonal.  Technically, this indicates singularity, but it is often
  49.      caused by improper arguments or improper setting of LDA .  It will not
  50.      occur if the subroutines are called correctly and if SBGCO has set RCOND
  51.      .GT. 0.0 or SGBFA has set INFO .EQ. 0 .  To compute  INVERSE(A) * C
  52.      where  C  is a matrix with  P  columns
  53.         CALL SBGCO(ABD,LDA,N,ML,MU,IPVT,RCOND,Z)
  54.         If (RCOND is too small) GO TO ...
  55.         DO 10 J = 1, P
  56.         CALL SGBSL(ABD,LDA,N,ML,MU,IPVT,C(1,J),0) 10 CONTINUE LINPACK.  This
  57.      version dated 08/14/78 .  Cleve Moler, University of New Mexico, Argonne
  58.      National Lab.  Subroutines and Functions BLAS SAXPY,SDOT Fortran MIN0
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.